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Detection of the quantum fluctuations by conventional methods meets certain obstacles, since it 
requires high frequency measurements. Moreover, quantum fluctuations are normally dominated by 
classical noise, and are usually further obstructed by various accompanying effects such as a detector 
backaction. In present work, we demonstrate that these difficulties can be bypassed by performing 
the cross-correlation measurements. We propose to use a pair of two-level detectors, weakly coupled 
to a collective mode of an electric circuit. Fluctuations of the current source accumulated in the 
collective mode induce stochastic transitions in the detectors. These transitions are then read off 
by quantum point contact (QPC) electrometers and translated into two telegraph processes in the 
QPC currents. Since both detectors interact with the same collective mode, this leads to a certain 
fraction of the correlated transitions. These correlated transitions are flngerprinted in the cross- 
correlations of the telegraph processes, which can be detected at zero frequency, i.e., with a long 
time measurements. Concerning the dependance of the cross-correlator on the detectors' energy 
splittings El and £2, the most interesting region is at the degeneracy points ei = ±£2, where it 
exhibits a sharp non-local resonance, that stems from higher order processes. We find that at 
certain conditions the main contribution to this resonance comes from the quantum noise. Namely, 
while the resonance line shape is weakly broadened by the classical noise, the height of the peak is 
directly proportional to the square of the quantum component of the noise spectral function. 

PACS numbers: 72.70.-|-m, 42.50.Lc, 74.50.-|-r, 73.23.Hk 



I. INTRODUCTION 

It has been recently understood that the noise phenom- 
ena are not necessarily something detrimental in physical 
experiments, but instead, they may carry a useful new 
information about the underlaying processes This in- 
formation may be difficult to extract from the measure- 
ments of average quantities, therefore properties of the 
fluctuations become themselves a valuable subject of re- 
search. Particularly, one may study higher irreducible 
moments (or cumulants) of physical quantities.'^ It is es- 
pecially interesting to detect essentially quantum prop- 
erties of noise, such as non-symmetrized correlators, re- 
flecting the quantum non-commutativity;^ 

To motivate this interest, it is useful to introduce an 
example of a quantum mechanical operator j{t), repre- 
senting a fluctuating electric current. Suppose this cur- 
rent is measured by a classical ammeter of certain band- 
width, and the raw measurement data is stored as a se- 
quence of numbers {j{tk)). The classical properties re- 
lated to the symmetrized correlators, such as the noise 
power S'sym(w) = J dTe'icp{iu!T){{5j{t),5j{t + t)}), or 
higher order symmetrized correlators, can be extracted 
from the raw data in postprocessing. By contrast, the 
quantum (anti-symmetrized) parts of the correlators can 
not be found in principle from postprocessing of such raw 
data, since this information is lost in the measurement of 
the average values {j{tk))- Instead, these quantum prop- 
erties may be inferred from the measurement of more 
complex quantities and systems, occurring naturally, or 
engineered on purpose. However, in this kind of mea- 
surements, it may not always be immediately clear which 
particular quantity is measured. Specifically, there is cer- 
tain ambiguity concerning the questions of the operator 



ordering in the complex quantities, and of the process of 
their reduction to the classical values. Thus, a careful 
approach requires the knowledge of the detailed model of 
the detector and of the process of measurement. 

Hence one arrives at the notion of a mesoscopic on- 
chip detector^Si^— It is a part of the measurement appara- 
tus, which interacts directly with the system, and is re- 
sponsible for transforming the quantum information into 
a classical signal. As an example, the two-level detec- 
tors have been studied theoretically and implemented in 
experiments,— and their operation is now quite well 
understood. Such a detector consists of a double quan- 
tum dot, or a similar structure, which has two energet- 
ically relevant quantum states (see Fig. [ij, that have 
different spatial charge distributions. Then a noninva- 
sive QPC electrometer—"— may be used to read out the 
state of the two-level detector, and the resulting signal 
can then be amplified by conventional means. In a prop- 
erly adjusted operating regime, the two-level detector is 
weakly coupled to the mesoscopic system, and at short 
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FIG. f : A schematic representation of the two-level detector 
as a double well structure with the level splitting e. Stochas- 
tic transitions in the detector are induced by the tunneling 
coupling A and assisted by photons coming from the collec- 
tive mode. The state of the detector is monitored by the QPC 
electrometer, located nearby one of the wells. 
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times they evolve together quantum mechanicahy. Due 
to the weakness of couphng, fluctuations in the system 
induce rare stochastic (nonadiabatic) transitions in the 
detector. Because of the noisy nature of these fluctu- 
ations, the state of the detector becomes decohered to 
a statistical mixture.— i^^""— In this case, the QPC elec- 
trometer effectively senses already classical state of the 
two-level detector, and thus, it actually satisfies our def- 
inition of an on-chip detector. 

While the on-chip detector approach should at least 
clarify as what exactly is measured, the extraction of the 
information about the quantum fluctuations may still be 
challenging for some general reasons. The quantum ef- 
fects often appear as small corrections to classical con- 
tributions, thus a high relative accuracy may be needed. 
Another source of complication is the fact that the system 
of interest is subjected to the perturbations induced by 
the measuring device itselff^ and other extrinsic sources 
of noise. Therefore, some advanced techniques have to be 
employed to carefully extract the useful information. One 
such technique, that demonstrated certain success at iso- 
lating the properties of the measured system is a so called 
cross-correlation techniqucj ^^'^" The main idea is that 
two (or probably more) detectors are used to measure 
certain fluctuating signal from the same system. Then, 
since the detectors are meant to be independent, any lo- 
cal processes at one of the detectors should not lead to 
cross-correlations in the fluctuations measured by differ- 
ent detectors. Therefore, the measurement of the cross- 
correlations of the detectors' outputs, gives certain level 
of protection from the local unwanted sources of noise, 
and from the detectors' backaction to the system, and 
thus may enhance the accuracy of the experiments.—"— 

In present work we consider using the cross-correlation 
technique to gain access to quantum fluctuations of cur- 
rent. We propose the following measurement setup: a 
mesoscopic system, source of current noise, is incorpo- 
rated into an electric circuit, so the fluctuations of its cur- 
rent are accumulated on a capacitor (see Fig. [5]). An elec- 
tric charge on the capacitor plays the role of a so-called 
collective mode, to which a pair of two-level detectors are 
coupled. We assume the weak coupling regime, when the 
evolution of the detectors' state can be described by the 
master equation for the average state occupations. Then 
fluctuations of the collective mode induce rare stochastic 



to detectors 







SOURCE 


I R 





FIG. 2; The equivalent scheme of the measurement electric 
circuit coupling the noise source j{t) to detectors via the col- 
lective mode Q. 



transitions in the detectors. These transitions generate 
two telegraph processes in the outputs of corresponding 
QPC electrometers, (see Fig. |4|) and the cross-correlator 
of these outputs has to be measured in the markovian 
(long time) limit. 

small A: perturbative I cumulant 

small a: expansion I resummation 
■ 

Tc < Td <. Ts Tm t 

■ I 

coherent evolution | | stochastic transitions 

FIG. 3: A diagram of the hierarchy of the time scales in 
the measurement process: the noise correlation time Tc, the 
decoherence time Td, switching time of the detectors Ts, and 
and the proposed measurement time Tm- 

Throughout the paper we rely on the concept of the 
time scale separation. In our model there are several 
independent small parameters (such as the coupling con- 
stant a, and the tunneling amplitude A) that produce the 
following grid of time scales (see Fig. [3]) . The smallest 
scale is the noise correlation time Tc , at which the coher- 
ent quantum-mechanical evolution of the joint system of 
the noise source, the circuit, and the detector takes place. 
However, starting from the decoherence time Td Tc, the 
noise source plays role of the heat bath for the detector. 
Then the evolution of the average occupations of the de- 
tector states is described by the master equation. These 
occupations vary on the characteristic time Ts ^ Td of 
the order of the inverse transition rates. Finally, the 
longest time scale Tm ^ Ts is the markovian limit of the 
telegraph processes, where the cross-correlator should be 
measured. 

To describe the quantum mechanical evolution of the 
joint system of the collective mode cj) and the two-level 
detectors on short time scales t < r^, it is rather more 
adequate to think of the two detectors in the space of four 
states 1 11), 1 12), |21), |22). A calculation of the reduced 
density matrix for the detectors involves averaging over 
the fluctuations of the collective mode (f>. To the nth 
order of the perturbation expansion with respect to tun- 
neling A, this requires finding averages of the Keldysh 
ordered products of the corresponding number of expo- 
nential phase operators e*"*^ (so called vertex operators) . 
In the weak coupling regime a <^ 1, we use the cumulant 
expansion, since every next cumulant enters with one ex- 
tra power of the small coupling constant, and we limit 
ourselves to the third cumulant. The main contribution 
comes from the times as long as the decoherence time 
Td, where the cumulants can be expressed in terms of 
the zero-frequency expansion of the noise spectral func- 
tions. Quantum corrections to the classical long time 
asymptotic of the cumulants are small, and thus, they 
can be taken into account perturbatively in the coupling 
constant a. 

On the time scales t > Td, we find the transition rates 
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perturbatively in A. We show that the most interesting 
effects do not appear on the level of the standard P{E) 
theory,—!^ which accounts tunneling to the lowest (sec- 
ond) order, only. Notice, that some of the transitions 
between these states, such as |11) ^ |22) or |12) ^ |21), 
that arise in the 4th order in A, directly correspond 
to simultaneous, correlated switching of the detectors 
and lead to cross-correlations at long times (see Fig. 
However, computation of the transition rates in the en- 
ergy representation encounters well known divergences 
in higher orders of the perturbation theory. We avoid 
those divergences by considering directly the time evolu- 
tion of the detector states. Namely, at times longer then 
the decoherence time t > Td, but much shorter than the 
switching time t <^ Ts, the master equation description 
suggests linear in time drift of the occupations probabil- 
ities from the initial distribution. On the other hand, in 
the fourth order we find quadratic in time terms. We 
identify them with the reducible contribution generated 
in the perturbation expansion and, accordingly, find the 
transition rates by extracting the irreducible part from 
the long time asymptotic of the occupation probabilities. 
As typical for the perturbation theory in higher orders, 
immense numbers of terms are generated. Nevertheless, 
we managed to find all our results analytically 4^ 

Note, that the transition rates may be measured 
from the time-resolved observation of the telegraph 
processesii^"— Then, by fitting their dependence on the 
controllable parameters of the system, such as energy 
splittings and coupling constants, one may try to infer 
some of the noise properties. The drawback of this ap- 
proach is that the large amounts of real-time data need 
to be recorded and analyzed. This also limits the possible 
measurement pace, since the real-time switching resolu- 
tion is required to extract the transition rates. Instead, 
we expect that by measuring directly the cross-correlator 
of the two telegraph processes on the time scale Tm ^ Ts 
longer than the switching time, one can considerably sim- 
plify the implementation of experiments. To evaluate 
the resulting cross-correlator of the telegraph processes, 
we generalize for the case of two detectors the approach 
that has been proposed in Refs. [4l| and to study 
the statistics of bistable systems. We present an exact 
general result, which is convenient to use if the transi- 
tion matrix is symmetric, i.e., for the classical noise. For 
the case of quantum noise we also develop a perturba- 
tive in tunneling calculation which is better suitable for 
analytical computations. 

Finally, we analyze the physically different contribu- 
tions to the cross-correlator from classical and quantum 
noise. Note, that although the two detectors are not di- 
rectly interacting one with each other, an effective non- 
local interaction emerges between the detectors due to 
their coupling to the collective mode [see Eq. dS])]. The 
effective interaction does not depend on the properties of 
the noise, but it leads to the trivial second order contri- 
bution to the cross-correlator. This contribution is not 
of much interest, since it can be activated by subjecting 
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FIG. 4: An example of the two telegraph processes gener- 
ated by (partially) correlated switching of the two detectors. 
One of the mechanisms of cross-correlation is illustrated by 
correlated transitions |11) <=^ |22) (red), and anti-correlated 
1 12) ^ 1 21) (blue), both governed by the anti-diagonal ele- 
ments of the transition matrix. The net cross-correlator may 
be estimated simply by counting the number of red lines, mi- 
nus the number of blue lines. 



both, or even any one of the detectors to local classical 
noise. However, since it is proportional to the strength 
of the effective interaction, it can be minimized by tun- 
ing the parameters of the circuit. By contrast, the more 
interesting kind of contributions arise as a consequence 
of the correlated response of the detectors to the fluctu- 
ations of the collective mode in the higher (fourth) order 
processes. The distinctive feature of these contributions 
is a sharp peak of the non-local resonance at the degen- 
eracy point, where the detectors' level splittings satisfy 
£i ± £2 = 0. We note, that the non-local resonance con- 
tains both classical and quantum parts of noise. Remark- 
ably though, the classical part in the 4th order is also pro- 
portional to the effective interaction strength and thus 
may be rendered small, while the purely quantum con- 
tribution to the cross-correlator survives if the effective 
interaction is "switched off". The corresponding condi- 
tions (see Sec. IVIC|) can be achieved, basically, if the 
noise temperature is sufficiently high, and the coupling 
constants are made small enough. 

The paper is organized as following; we present the 
model of two-level detectors in Sec. [Til Then, in Sec. IIIII 
we give an overview of the counting statistics approach, 
generalize it for the cross-correlations, and present two 
methods of evaluating the cross-correlator, usable for 
classical and quantum noise respectively. In Sec. IIVI and 
|V]we show how the cumulant expansion is applied to cal- 
culate the averages needed to find the transition rates. 
Finally, in Sec. IVII we present the results for the classi- 
cal and quantum noise cases, and analyze the conditions 
needed to access the required measurement regime. 



II. TWO-LEVEL DETECTORS 

An electric circuit may be modeled'^^ by a set of 
bosonic fields, 4>k, and their conjugated "charges" qk- We 
wish to single out one of these fields, a so called collective 
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mode Q and its phase 0, which are linearly coupled to 
the rest of the fields and enter quadratically to the corre- 
sponding Hamiltonian of the circuit. One can show that 
any such Hamiltonian may be transformed to the form, 
where all the couplings are carried by the phase (j) only: 

'He ^ +'Hn{<P,qk,(l)k)- (1) 

This model is sufficiently general to describe any cir- 
cuit with the usual linear elements, such as resistors, 
capacitors and inductances, and can also include such 
mesoscopic elements as tunnel junctions by adding non- 
quadratic potentials to the Hamiltonian. 

Let us assume that the collective mode Q is linearly 
coupled, with a dimensionless coupling strength a, to a 
two-level detector represented by the Hamiltonian 

•Ho = ^cTz + Aa^ + ^cTzQ, (2) 

where e is the energy level splitting, A ^ e is a weak 
level mixing and we use the units e = h = I through- 
out the paper. The constant A can be also understood 
as the tunneling amplitude between the levels, or as the 
quantum level broadening. This type of quantum detec- 
tors have been considered in a number of works and have 
been implemented experimentallyj^^— 

In what follows, we treat the term Aa^ as a smallest 
perturbation. We, therefore, render the total Hamilto- 
nian = T-lo + "He in a more convenient form by per- 
forming the following transformation: 

H' = e*^''="*-He~*^'"-'^. (3) 

This transformation affects only those terms that do not 
commute with CTz or 0. Since [(j), Q] = i, it shifts the 
charge Q by —aa^/^, canceling the linear coupling term 
aOzQI^C and bringing interactions in the form of oper- 
ators e^'""^, 

where the energy is also shifted by a constant —o?I^C. 
We switch to the interaction picture with the time de- 
pendent tunneling Hamiltonian 

/ gi[</'(i)+et] \ 



Note that if i?c = the tunneling Hamiltonian (|TT|) re- 



This suggests for a perturbative expansion in powers of 
A. To justify this expansion, we assume that the A has to 
be the smallest energy scale in the system. Particularly, 
it must be smaller than the level broadening introduced 
by noise. 

A. Pair of two-level detectors 

Now, let us consider two such detectors, both coupled 
to the same circuit via the charge Q. 

j=l,2 j = l,2 

(6) 

where we denote quantities belonging to each one of the 
detectors with an additional index j = 1,2. In this case, 
the analogous transformation, 

leads to 

+ He+ l^aWai^). (8) 

Thus, besides changing the energy by a constant, it also 
generates the cross-term with 

Ec = aia2/2C. (9) 

The cross-term is nothing but the Coulomb charging en- 
ergy of the capacitor for the screened detector charges, 
and represents the effective non-local interaction between 
the detectors mediated by the circuit collective mode. At 
the degeneracy points Si = ±62 this interaction leads to 
the quantum avoided-crossing level splitting of a value 

Ae = 4AiA2E^/el (10) 

Finally, after switching to the interaction picture, the 
tunneling Hamiltonian for two detectors takes the form 



(11) 



duces to a tensor sum Hi = U^P «) E + E (g) 'Hf\ We 



Hi 



Q ^^^-iax4>-i{e^+Ec)t ^^g-ia24>-i{e2+E^)t q 

I 
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will see that the presence of Ec leads to a trivial mech- 
anism of cross-correlations even in the presence of only 
local noise. However these cross-correlations vanish with 
small Ec, and thus may become dominated by some more 
interesting phenomena, such as effects of quantum fluc- 
tuations in higher order processes. 



B. Time evolution 

To study cross correlations in the detectors' output, 
we wish to consider the long-time limit, where evolution 
of the detectors' states can be described by the master 
equation 



MP 



(12) 



for the occupation probabilities P = (pi, . . . ,^4). Such 
description is valid if the decoherence time — 
{a^R^S)~^ is much shorter, <C A~^, than the quan- 
tum time scale associated with the quantum level repul- 
sion of a single detector, or Td <C (Ae)^^ in the degen- 
eracy point. This conditions are equivalent to a require- 
ment that the classical level broadening due to the noise 
is always stronger then the quantum avoided-crossing 
level repulsion. Or, in other words, one may say that 
the stochastic switching time Ts = l/(A^Trf) is longer 
than the dephasing time, 3> t^. 

The time dependence of the reduced density matrix of 
the detectors is 



where 



p(i)-Tr,[C/p(0)C/t], 



U^Texp{-i J dt'Uiit') 



(13) 



(14) 



is the interaction picture evolution operator, the ini- 
tial condition is represented by the full density ma- 
trix of the system of detectors together with the circuit 
/5(0) = p(0) X Pc(0), and Trc means averaging over the 
circuit degrees of freedom. 

As soon as p{t) is found, one observes that the off- 
diagonal elements of p decay exponentially over the time 
Td- Therefore we can concentrate our attention on the 
probabilities 

Pfc(i)=PfefeW = Tr[p(0)xp,(0) U^\k){k\u\, (15) 

where the trace is over all the degrees of freedom of the 
system. The equation (jlSp may be recast in the form 



P(t) = [E + ^(t)]P(0), 



(16) 



where the elements of the time dependent transition 
probability matrix J^{i) read 



= Tr, 



pM {k\U^\l){l\U\k) 



(17) 



They may be found from (|15p as the perturbative expan- 
sion M.{t) = Yln-^riit) with respect to the tunneling 
Hamiltonian pT|) . Comparing the expression (|16l) with a 
solution of the equation dH]), P{t) = exp{Mi}P(0), one 
finds that the master equation matrix can be obtained 
from the long time asymptotics of the irreducible part of 
M{t) 



M = lim - log [E + Mit)]. 

t^oc t 



(18) 








1 



t 



FIG. 5: Typical, quadratic in time dependence (red line) of 
the elements of bare A^4(t), and linear in time dependence 
(blue line) of corrected jCi4{t) — Ai2(t)/2. Time is indicated 
in units of the decoherence time Td. 

Note, that due to the specific structure of (dJ), only 
even powers of Aj are present in the diagonal elements 
of the perturbation series for the density matrix p^ . 
Therefore, to the fourth order in Aj, one can write 
M{t) = M2it) + M4{t). The time dependence of an 
arbitrary element of the M2{t) appears to be quite sim- 
ple: 



(19) 



where Ci is a constant, and the first term Co{t) is a de- 
caying oscillating function of time. This decaying term 
only represents the artefact of the chosen decomposition 
of the initial condition p = p x pc, which, due to the 
detector-circuit interaction, does not describe a station- 
ary state. However, the stationary, linear in time behav- 
ior is restored after the decoherence time Td , as illustrated 
by the decaying oscillations in the Fi^. [S] 

The elements of the next term, Ai^^t), may have a 
more complex structure: 

{Miit)}^^ = C'„it) + [C[ + D[it)]t + C'^t\ (20) 

where the functions C'^it) and D[{t) are analogous to 
Co{t), the constant C[ is analogous to Ci, and there is a 
quadratic in time term G^t^ . This last term is a reducible 
part of the A^4(<), which is exactly eliminated by the 
logarithm in (IT51) . Particularly, to the 4th order one has: 



M 



lim — ^ 

t->oo t 



M2it) + Miit) - -Mlit) 



(21) 



This expression provides a conclusive point in the calcu- 
lation of the transition matrix M from the explicit time 
dependence of the occupations P{t) . 
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III. FULL COUNTING STATISTICS OF THE 
QPC CURRENTS 

In this section we present the method for calculating 
the statistics of a QPC current using the generating func- 
tions approach. We first recall the simpler case of a single 
detector— and later generalize it to the case of two de- 
tectors. Consider a QPC as a charge detector tuned to 
detect the state of the two- level system (see Fig.[T]). The 
statistics of the current passing through the correspond- 
ing conductance levels A: = 1,2 of the QPC is described 
at long times by the moment generating function 



Hfc(A)i 



(22) 



where fk (n) is the probability that n electrons are trans- 
ferred while the levels k is occupied, t is the total time of 
the measurement, and -fffe(A) is the function generating 
current cumulants: 



{{in = 



One can see that gk{X,t) satisfy the equations 
5fe(A,t) = i^fe(A)gfc(A,^)■ 



(23) 



(24) 



This equations may be modified in order to take into ac- 
count mixing of the current channels induced by switch- 
ing of the detector. Introducing G(A,i) = (51,52), the 
extended master equation reads: 



G(A,i) = WG{X,t), 



(25) 



where the matrix W = H + M is a. sum of the cumulant 
generating functions matrix 



H 



ffi(A) 
i/2(A) / ' 



and the transition matrix 
M = 



(26) 



(27) 



where the transition rates are denoted as r±. In the long 
time limit, the statistics of the QPC current is given by 
the generating function 



H{X)^ \mi ilogV5fc(A,t), 

t^OC t ^ — ^ 



which is the largest eigenvahie of the matrix 



H{X) = H,-r, + V(i?d + rd)2 + r_r+, 



(28) 



(29) 



where H^^a ^^{Hi±H2), and Ts,d = i(r+ ±r_). Cur- 
rent cumulants are then found as derivatives of this gen- 
erating function 



(/")) = 



aA" 



(30) 



Taking into account that 



dH{X) 



dX 



A^O 



E 



dH dHi. 



dm ax 



A^O 



dH 



dm 



(31) 



we obtain the expected result for the average current: 

{I) = ^-^ + ^^Mh)+P.{h), (32) 
and, similarly, for the noise 

= {h)y+pi{{in (33) 

Here (Ik) is the average current at the conduction level 
k of the QPC, p2,i — r±/(r+ + r_) are the stationary 
level occupations, and ((/|)) is the corresponding zero 
frequency noise power of the QPC. 



A. Cross-correlator 

The above scheme can be easily generalized to a pair of 
two such detectors. Without level mixing, the joint prob- 
ability distribution for the numbers of electrons passed 
through each QPC factorizes: fki{ni,n2) = fk{ni)fi{n2), 
therefore the corresponding generating functions acquire 
the form 



gkiiX,T^,t)^gi'\X,t)gP{7it), 



(34) 



with different counting variables A and rj for the two 
QPCs respectively, and satisfies the equation 



5fci(A,77,0 = [Hj^'{X)+Hl^'{7j^gki{X,V,t), (35) 

where the upper index (j) indicates the properties of the 
detector j. 

Switching of the detectors is accounted for in the ex- 
tended master equation, G — WG, with the matrix 
that has the same structure as before: W = M + H, 
where M is the transition matrix. In the basis of states 
(|11), |12), |21), |22)), the diagonal matrix H of the cu- 
mulant generating functions is: 



(2), 



H ^ m 



E + E(g)i?2, 



where Hj are the single detector matrices ([26] 



H 







(36) 



(37) 



as functions of corresponding counting variables A or rj. 

As before, current cumulants may be found as deriva- 
tives of the largest eigenvalue H{X,r]) of the matrix W. 
Particularly, the current cross-correlator reads 



{{hl2)) = d^d„H{X,f^)\^^^^^. 



(38) 
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Since there is only one derivative with respect to each 
counting variable, the cumulant generators Hj: can 
be replaced by their corresponding average currents 
ijf )(A) w A(4^^) and iff ^(r?) « 77(/f Thus the prob- 
lem of finding the cross-correlator of QPC currents can be 
reformulated as finding a second order correction to the 
largest eigenvalue of the matrix M from the perturbation 
H. 

Note, that due to the conservation of the total proba- 
bility Mjk = 0, the transition matrix is degenerate, 
i.e., one of it's eigenvalues mo — 0, and all the others are 
negative, rrii^o < 0. Therefore we can use the general re- 
sult of the perturbation theory that i/(A, 77) can be found 
as a second order correction to the largest eigenvalue mo : 



mo - rrii 



{mo\HM*H\Tno), (39) 



where M* = '^i^Qi^/'nii)\mi){'mi\ is also known as a 
Moore- Penrose pseudoinverse of the matrix M; jm^) are 
the eigenvectors''^ of M. 

To find the cross-correlator we need to calculate the 
derivative in Eq. psp . Since the counting variables only 
enter the matrix H, it is convenient to introduce the fol- 
lowing operator 



dxdri 



(H\mo){mo\H) 



(40) 



so that the cross-correlator may be expressed as follow- 



ing: 



((/1/2)) = -Tt(V-M*). 



(41) 



It appears also, that the cross-correlator is always 
proportional to the differences of the current levels 



r(2)\ 



(I- 



(2)\ 



therefore 



A/i = {r^'>) - {I\'>) and A/2 = (/; 
it is natural to define a normalized cross-correlator, that 
has the dimensionality of time. 



X 



{{hh)) 
A/1A/2 



Tr(X' • M*) 
A/1A/2 ■ 



(42) 



Finally, we note that in practice, for systems with suf- 
ficiently general transition matrix M, the approach de- 
scribed above typically leads to cumbersome expressions, 
in particular, for the eigenvalues and the eigenvectors of 
M . Therefore, we present a different way to evaluate the 
cross-correlator for a particular class of transition matri- 
ces. 



B. Perturbative approach 



neglected, iJc — >■ 0. In this case one can easily find the 
cross-correlator perturbatively in tunneling amplitude A. 
The transition matrix can be represented as 



M = Afo + 5M, 



(43) 



where Mo is evaluated to the 2nd order in tunneling and 
represents completely uncorrelated detectors, while 5M 
is a small correction of the 4th order. The matrix Mo = 
Ml (g) E -I- E (g) M2 is a tensor sum of the single-detector 
transition matrices Mi and M2 j defined according to Eq. 
dSZl) as 



Mj = 



-r 



r 



(j) 



-r 



(j) 



(44) 



The extended master equation matrix then reads 

W ^M + H = Wq + SM, (45) 

where, the matrix Wq — Mq + H retains the same tensor 
structure: Wq = W^i(g)E-|-E(8)Vi/2, and Wi and W2 are the 
single detector matrices as in ([^5)) . The matrix Wq can be 
diagonalized explicitly, and it's largest eigenvalue is, ob- 
viously, the sum of the corresponding eigenvalues (|29l) for 
each detector separately, -ffo(A,?y) = iJ^^^(A) + H^^^r]). 
Thus, as expected for the uncorrelated detectors. 



(46) 



where the index stands for the lowest order contribu- 
tion. We see that any cross-correlations may only be 
generated by the correction matrix SAI, which brakes the 
tensor sum structure of W .'^^ 

As before, we need to find the largest eigenvalue 
H{\,rj) = Fo(A,77) +SH{X,r]) of W, where SH{X,r]) is 
a correction due to the perturbation 5M. Note, that in 
contrast to the general case, the counting variables are 
not contained in the perturbation SM, but are absorbed 
into the matrix Wq (which makes Wq non-degenerate). 
Therefore, it is sufficient to use the first order perturba- 
tion theory: 



(47) 



where |mo) is the eigenvector of Wq corresponding to the 
eigenvalue HqIX^tj). 

Finally, the cross-correlator is given by 



iihh)) ^ dxd^dHiX,ri)\^^^^^ 
and analogously to (pij) . we find that 
d^SH 



(48) 



{{hh)) = E 



(2) 



(49) 



We wish to consider a particular case, when the ef- 
fective interaction Ec in the Hamiltonian (|lip can be 



Although the full expression for SH is somewhat cumber- 
some, it has an important property that the derivatives 
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with respect to h[^^ and Hi^\ all proportional to the 
same quantity Xp, 

d^SH 



X.,7J=0 



= (-1)'=+%. 



(50) 



Therefore, the cross-correlator ([^^ acquires the same 
form as in (|42]) and can be found using one of the re- 



lations (l50l 



2p(l)p(l)p(2)p(2) 

(r(i) + r(i))3(r(_2) + rf )3 



^<5M,,%, (51) 



where to present a compact expression for Xp, we have 
introduced the matrix of the coefficients K,-,-: 



J 



K 



n(2) 



r 

Ar<i)Ar(^: 



2r\ 



T-(l) 



- 1 



Ar<^>Ar<^' 



2r^ 



+ 

-1 

p(2) 



' + 



(52) 



and the notation AF^^' = v'f^ - tL^^ is used, with t'£ 
being the transition rates to the second order of the per- 
turbation theory. 



IV. CUMULANT EXPANSION 



To find the transition rates p3|) . we expand the transi- 
tion probabilities (jl7p in the perturbation series with re- 
spect to the tunneling Hamiltonian (fTT|) . We will hence- 
forth use a simple angled brackets to denote the weighted 
average (...) = Trc[pc---], which is used in ([T7|. By 
substituting the evolution operator (fT4l) . expanded per- 
turbatively with respect to the tunneling Hamiltonian 
dTTI) , into the expression (IT7)) for the transition probabil- 
ity matrix, one can show that the elements of Ai (t) may 
be expressed in terms of time integrals of the Keldysh 
ordered correlation functions of the following type: 



(53) 



where = 4'{tj) for simplicity, the 7j take values of ±ai 
or ±a2, and the generating function 



k 



■'/'fe)) (54) 



is an expansions in terms of the Keldysh ordered cumu- 
lants Tk((^!>i . ■ • 4>k)} of the order k. 

We do not need to know all the correlators that appear 
in the expansion of the density matrix, but only those 
that contribute to the master equation through the ma- 
trices M2{t) and M4{t) from (|2T|) . In the second order 
term A42{t), due to the structure of the tunneling Hamil- 
tonian (jlip . we find only the terms proportional to 



or A2, thus "originating" separately from each detector. 
In this case, we need the correlator of the form 



(55) 



('2) 

where the second and third cumulants, J2 (^1,^2) ^nd 
--^3 (^1,^2), are found by expanding both sides of the 
equation (|53p in 

2 

4^\tiM)^\{4>\-'24>i4>2 + (t>t), (56) 

3 

4\ti^t2) ^\{<t>l- i<t>l4>2 + iMl - 4>f)- (57) 

By contrast, in the fourth order term Mi{t), there is 
more diversity. First of all, there are terms proportional 
to A^Aj. Such terms are of our main interest, since they 
represent the non-local processes, involving both detec- 
tors. Thus, we are interested in correlators 



0*7101 i7404\ _ -J^,^''(tl;---,ti)-iJ^^''{ti....,ti) 



(58) 



where the four parameters 7j take four different values 
of ±ai and ±a2. Again, expanding both sides of the 
equation ([55)1 . one finds 



k k<l 



and, in contrast with the result (|57p . the third order cu- 
mulant 

i i<j 

+ 6^ ^ajjk4'i4>j(t)k'^ (60) 



i<j<k 
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contains the three-point correlators 
three different times. 



taken at Taking the integrals over w, we find the long time 

asymptotic behavior of the correlator (pS)) as 



V. CURRENT AND PHASE CORRELATION 
FUNCTIONS 

From the equation of motion cj) — Q/C generated by 
the Hamiltonian ([T]), it follows that the phase is an 
integral of the voltage on the capacitor. Since the con- 
stant component of voltage can be absorbed into the level 
splitting energy e, we include only the fluctuations 5V , 



4){t) = / dt'5V{t'), 



(61) 



where instead of integrating over time from minus in- 
flnity, we have introduced the regularization with some 
finite time — T, which is larger than any other timescale 
in the system. 

The voltage fiuctuations are related to the current fiuc- 
tuations by the circuit impedance Z{ui)^ 



5V{u:) = Z{uj)3{u). 



(62) 



Combining the last two equations together, one obtains 

(63) 



J 2tt (-iw) 

where, due to the regularization, a pole at w = is can- 
celed by a zero in the numerator. We further assume the 
limit of fast circuitp^ i.e., the main contribution to the 
phase correlators comes from the times much longer than 
the circuit response time trc = _RC, and the impedance 
may be considered a constant Z{lu) — R. 

We define the current spectral functions S2{ijj) and 
Sz{loi^uj2) through a non-symmetrized two- and three- 
point current correlation functions 



(j(wi)j(w2)) = 2tt6{uJi +UJ2)S2{iOi), 



(64) 



{i{^i)i{^2)3{<^3.)) = 27r^(wi +UJ2 +uj'i)S'i{uJi,uj2), (65) 

where the delta-functions are the manifestation of the 
time translation invariance, since we consider a station- 
ary process. Using these definitions, the phase correlators 
are expressed as following: 



b{yi\ -t- W2) 



S2{^i) n 



(66) 



and analogously. 



(27r)3 



X 5'3(a;i,W2) J]^ 

i=1..3 



5{uJi + U]2+ 1^3) 



- iS' sign(ii - tj) . (68) 

Here only the zero-frequency components of the spectral 
function contribute to the correlator: 



(69) 



where S = 6*2(0) is the classical noise power, and 
S' = dS2{uj) / duj\i^^Q is the quantum part, representing 
non-commutativity of the current operator. 

The correlator (j68p has a divergent term proportional 
to T, which, however, does not contribute to any phys- 
ically meaningful quantity. For instance, the cumulants 
([55| and read: 

J2{tl,t2) = ^2~(^'^' " " ' ^'^'"^^^ " ^'^)' ^'^^^ 

while the cumulant from the fourth order can be rewrit- 
ten in the similar form: 



J. 



(4) 



ih. 
i?2 



,^4) 



^(^S\U -tj\- iS' sign(<, - t^)). (71) 



A weak coupling regime, considered throughout the pa- 
per, implies that either the detector to circuit coupling 
is weak (small aj ^ 1) or the impedance R is small. 
Then the entire J2 is small by the factor Uia^R^, and 
the contribution to the time integrals of correlators ([55)) 
and ([58l) that enter the perturbation expansion for M.{t), 
see Eq. (fT7|) . comes from large differences of times of or- 
der \ti — tj\ ^ [oiiOijR? S)~^ . This justifies the use of 



asymptotic expressions of the type ([70)1 and ([71]) and 
clarifies the assumption of the fast circuit trc ^ ''d = 
{a^R^Sy^. 

Analogously, for the three point correlator (|67|) . we 
expand the spectral function around zero frequencies. 



S'3(wi, W2) = 5'3 + S^^^L^x 



(72) 



and obtain 



{cl){h)<P{t2mh)} - ~R'S3[T + minitj)] 



iR' J2 S: 



;mm(<,)-i,)-- 



(73) 



However, the parameters S^^^^ and S^'^' are not inde- 
pendent. From the definition (j65p one can find that 
(67) (5^^^)* = = 2ReSP, therefore they are parameter- 
ized by two real numbers: the real and imaginary parts 



:(2) 
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(2) 

gain, although this correlators diverge with T, 
the third cumulant, such as (j57p . is regular: 



•^3^^ (^1:^2) 



6 



'+zlm5f^sign(ti-i2)]}. (74) 



Note, that in contrast to the second cumulant ([TO]), which 
depends on the absolute value |ii — ^2!, the classical part 
of the third cumulant depends on the difference (^i — 12), 
and because of this, it enters the transition rates the same 
way as the energies ej . This leads to a specific shift of the 
energies, discussed later. A similar expression is obtained 
for the third cumulant (|60l) to the fourth order, but it is 
not displayed here, because it is somewhat cumbersome. 



found explicitly. Note, that since the matrix M is sym- 
metric, it is parametrized by 6 elements of its upper-right 
triangle. Then, finally, the classical contribution to the 
cross-correlator is expressed in terms of those elements: 



32E2 



(mi2- m34)(TOi3- TO24) - (mi4- m23)J2i (^5) 



where J2i = "^12 + "n^u + "^'24 + "134 , and the sum in 
the denominator ^2 — ''^ij^kifnnm is over products 
of all combinations of three different elements from the 6 
independent elements of M. 



VI. THE CROSS-CORRELATOR: RESULTS 

Using the cumulant expansions such as ([701 . ([TT]) . and 
(ffi)) . we can find the correlators (I55|) and ([55|) . which we 
need to calculate the transition probability matrix (|17p. 
From the latter, we extract the master equation matrix 
(PT|) and, finally, find the cross-correlator using one of the 
results (|^^ or ([?T|) . We assume that the cross-correlator 
is a sufficiently smooth function of the parameters of the 
system, and thus, consider separately different contribu- 
tions, each being suppressed by some small parameter. 
This allows us to estimate ratios of these contributions 
and determine regions of the parameters where one or 
another contribution is dominating. We start with the 
case when the effective interaction Ec in the Hamilto- 
nian is neglected, and the noise is classical. Then the 
detector outputs become completely uncorrelated. We 
consecutively "turn on" the effective interaction and the 
quantum parts of the noise cumulants to investigate sep- 
arately each contribution. 



0.5 £2 




0.5 £1 



FIG. 6: The classical contribution Xd — X^j -I- X^^' to 
the cross-correlator in units of the dephasing time l/fi = 
Rns)-^ as a function of the detector level splittings 
El and £2 in units of the effective temperature T^h, at Aj = 

0.070, aj = 0.55 ± 0.05, R = 0.12, and C = O.S/fi. Note that 

(2) 

the second order part X^j does not exhibit any singularities 
at the non-local resonances ei ± £2 — >■ 0, and the sharp peaks 
come from the fourth order part X^j 



(4) 



A. Classical noise 

For the classical noise, the current operator j{t) may 
be replaced by a classical variable, which commutes with 
itself in different moments of time. Therefore, the spec- 
tral function 5*2(0;) becomes a symmetric function, and 
it's derivative at zero frequency is S' — 0. In the same 
way, the three-point spectral function (1551) has zero slope 

) = S'f ) = 0, and we set 53(^1,^2) = 5*3- 
The transition matrix M for the classical noise is sym- 
metric. Taking into account that the sum over columns 
of the elements of M is zero, one can show that the sta- 
tionary solution of the master equation, corresponding to 
the largest eigenvalue, is |mo) = ;|(1, 1,1,1) with equal 
occupations oipj = 1/4 for all the four states, which can 
be regarded as a limit of high temperature. The matrix 
(PO)) in this case is a constant independent on the param- 
eters of the system, and the pseudoinverse M* can be 



The expression (j75p gives certain interpretation of how 
different combinations of the transition rates contribute 
to the cross-correlator, particularly, if the transition rates 
in one detector depend on the current state of the other 
detector. For example, if mi2 > 77134, 77113 > 77124, 
then the positive cross-correlations are induced. One 
important observation is that if the effective interaction 
term in the Hamiltonian is neglected, then the transi- 
tion matrix acquires additional symmetry 777^^ = mki for 
i + l = j + k = 5 (symmetry relative to the anti-diagonal), 
and 77714 = 77123 . Then the cross-correlator is identically 
equal to zero. In other words, the cross-correlator for 
classical noise vanishes together with the effective inter- 
action strength Ec- 

If the transition rates are calculated only to the sec- 
ond order in Aj, the anti-diagonal elements of M are 
zero, and the result (|75p simplifies further. Substitut- 
ing corresponding transition rates, we obtain the classical 
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contribution, shown in the Fig. [6] 



B. Quantum noise 



where we introduce the classical level broadening, 



(76) 



(77) 



and recall that Ec is the effective non-local interaction en- 
ergy The typical maximum magnitude of this quan- 
tity, that is acquired at £i/£2 ~ ±aiAi/a2A2, can be 
estimated as 



X, 



(2) 



± 



E^„ 



±- 



(78) 



It contains competing small parameters: Aj that are con- 
sidered the smallest energy scales in the system, and ef- 
fective interaction strength Ec- It is especially remark- 
able that Ec is a property of the circuit itself and does 
not depend on the source of noise. In fact, the quan- 

(2) 

tity reflects rather the effect of the effective non- 

local coupling, than of the common noise source. If we 
replace the common noise source with two independent 
local sources, or even only one local source of noise, the 
cross correlations will persist. For the same reason, this 
classical contribution ((75)) is a well bounded function of 
the energies £i and 62, and it does not exhibit singulari- 
ties at the degeneracy points ei = ±£2. 

By contrast, if the transition rates are calculated up 
to the fourth order in Aj, then we obtain a correction 

xj^f\ which manifests resonance peaks at £2 = ±£1, that 
become more pronounced as the coupling constants ap- 
proach each other. Naturally, this is the most interesting 
region of parameters. Using notations ai^2 — azL da/2, 
£ = (£1-1- £2)72, and rii,2 ~ = we present 

the asymptotic form of the fourth order classical contri- 
bution to the cross-correlator 



X 



(4) 



ni2El 



im-^[{e2-eiY + nl^y 



where 



^^12 = -{ai~a2fR^S 



(79) 



(80) 



is the classical broadening of the non-local resonance. 
Here we assumed that Q.12 is larger than the avoided- 
crossing energy splitting A£ and the classical level broad- 
enings due to the local sources of noise. The height of 
the peak is equal to 



maxX 



(4) 
;1 



E'i 



at £1 



(81) 



In this section we concentrate on the effects of quantum 
noise, i.e., we now take into account the antisymmetric 
parts of the spectral functions 5*2 (w) and S'3(u;i,a;2). At 
the same time, we set Ec = 0, so the tunneling Hamil- 
tonian pT|) becomes a tensor sum, as well does the sec- 
ond order part of the transition matrix, Af2, but not the 
fourth order correction M4. This allows us to obtain an- 
alytic results using the approach described in the section 
piIBp . where we associate M2 with Mq in (|43)) . and the 
fourth order part M4 with the perturbation 5M . 



0.5 £2 




FIG. 7: Cross-correlations due to the quantum Gaussian com- 
ponent of noise over the classical background displayed in Fig. 
[6l in the same units, as a function of the energies ei and ei at 
01,0:2 > 0, develop a sharp resonance peak on the diagonal, 
the height of which is proportional to the energy squared. 



In the limit of weak coupling, the main contribution 
from the correlator to the transition rates comes 
from long times. This means that the first term in the 
cumulant, S\ti — t2\, is much larger then the second term 
iS"sign(ii — 12). This allows to include the quantum cor- 
rections perturbatively. Since the third cumulant con- 
tains one more coupling constant then the second cu- 
mulant, we expect the main contribution to come from 
Gaussian fluctuations. Surprisingly, however, the first or- 
der contribution of S' vanishes, and we have to keep the 
expansion at least up to S'^ . The fact that the Gaussian 
fluctuations do not contribute to the cross-correlator in 
the lowest possible order with respect to coupling con- 
stant means that we have to check if the contribution of 
third order correlations may enter with a lower power of 
the coupling constant. Unfortunately, this is not the case, 
and non- Gaussian correlations only make a small correc- 
tion to the Gaussian contribution, as both of them enter 
quadratically into the cross-correlator (|82|) . As a function 
of energy splittings, Xq exhibits the non-local resonances 
at £1 = ±£2. The exact expression is somewhat cumber- 
some, thus we do not display it, but concentrate on one 
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of the resonances and find the asymptotic form: 



645*2 



16(5')' + [4a' + 15((5a)2]i?2(Re5; 



(2)n2 



(82) 



Note, that in contrast to Eg . ([ST]) the hight of the reso- 
nance in the present case 



maxXn 



41712^2 



(83) 



is proportional to the square of tlie energy e, as can be 
seen in Fig. [T] 

To complete this subsection, let us emphasize, that the 
resonant quantum contribution to the cross-correlator is 
an essentially non-local effect. In the limit when the ef- 
fective interaction Ec is rendered vanishingly small, the 
quantum contribution becomes the only one present, that 
carries the manifestation of the non-local quantum cor- 
relation. 



Analysis of physical conditions for the quantum 
noise detection. 



One may notice that in the results presented above 
in this section we have not included the classical non- 
Gaussian noise 5*3. This is because the effect of this 
component of noise has been discussed earlier— and this 
effect is easy to describe: it consists solely in a uniform 
shift of the energies ej — ij = Sj + ajR^S^/G, and, plus 
to that, an additional small shift of 3aia2SaR^S3, ap- 
plies to the energy position of the non-local resonance 
£2 — £i- Obviously, this shifts become irrelevant if we are 
interested in the regions of larger energies Sj ^ flj . 




FIG. 8: The cross-section of the resonant peak of the quantum 
contribution to the cross-correlator as a function of Sa = 02 — 
ai and the energy difference £2 — ei illustrates how the peak 
becomes thinner and, finally, disappears as 5a — >■ 0. The gap 
width is determined by the local source noise, that also limits 
the height of the resonance. Same units as in previous figures 
are used. 



Another effect that we have studied, but that was, for 
the sake of simplicity, omitted until now, is the influence 
of the local noise sources. These local sources may be 
modeled by an additional fluctuating phases /3j'/>j added 
to the collective mode phase ajcj), including respective 
coupling constants l3j. Introducing corresponding noise 

powers one can find, that at the local resonances, 

these will simply add up to the classical level broadening 
Q,j ilj + 0jS\l].. More important, however, is their 
contribution to the classical level broadening of the non- 
local resonance: 0,12 — ^ P'^S\[l.. Thus, if we wish 
that this does not significantly wash out the resonance, 
then, besides of requirement of weakness of the local noise 

f^'j^iol ^ ' difference Sa = 02 — ai must be 

kept large enough. The effect of local noise is best seen 
in the Fig. |S]as a finite width gap between the two peaks 
at 5a > and at Sa < 0. Certain optimal value of the 
coupling constants may be chosen, that maximizes the 
magnitude of the non-local resonance. 

We have to determine if there is a range of physical 
parameters, where the quantum contribution (j83p to the 
cross-correlator is greater then the classical resonance 
contribution ([5T|) and the classical background (|78l) . In 
order to do this, we concentrate again on the non-local 
resonance Sj = e where the quantum effects are maximal. 
We recall the expression for the effective non-local inter- 
action strength ([9]), and our assumptions, made through- 
out the paper: (i) the level mixing constant is smaller 
than the the classical level broadening A <C i^j; (ii) the 
circuit is fast trc ^ ^^'^ (iii) the zero- frequency 

expansion for the spectral functions is valid for energies 
smaller than the noise temperature e < Teg. Demanding 
that the quantum contribution is larger then the classi- 
cal background, Xq > X^^ , and using (iii), we obtain 
the inequalities 



SaEc^ 

Toff > e > i~Tcff, 

a A 



(84) 



where we have substituted S/S' = 2Tcff. These inequali- 
ties are consistent when Ec/A < 1 or 



a^R < Atrc < 1, 



(85) 



i.e., if the coupling is made sufficiently weak. However, 
since the background is smooth, one can demand that 

the quantum contribution is larger than only the reso- 

(2) 

nant part of the classical contribution, Xg > X^j ' . This 
requirement leads to the condition 



Ec 



(86) 



which is weaker then (|85|). since A Q. This condition 
may also be reformulated as 

TcffTRc > 1. (87) 

Hence one can see, that although trc ^ flj^ = 
l/{a'^RTcs), the above inequality can be easily satisfied, 
in the weak coupling limit 1/a^R ^ 1. 
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Finally, we can also compare the cross-correlator with 
one of the auto-correlators Aj — ((/^))/(A/j)^, which 
can be found from ()33p and estimated as inverse of 
the typical transition rate, far from the local reso- 
nance Aj « e'j /{A'jU.j). Thus, even in the resonance, 
the cross-correlator to auto-correlator ratio, Xq/Aj ^ 
a^B? / (<5a)^ <C 1, is still much smaller then one. In other 
words, correlated transitions make only a small fraction 
of all transitions. 



VII. SUMMARY AND CONCLUSIONS 

Motivated by the challenges on the way of measur- 
ing the properties of quantum fluctuations, we have pre- 
sented in detail one possible approach to the problem, 
based on the cross-correlation technique. We use a con- 
cept of the collective mode, through which fluctuations 
of the current of the noise source are transferred to the 
on-chip detector. The latter consists of a pair of the 
two-level systems, monitored by QPC charge detectors. 
Thus, noise induced transitions generate two telegraph 
processes in the QPC currents, cross-correlator of which 
is supposed to be measured. The cross-correlation tech- 
nique allows one to avoid the effects of local noise sources 
and can provide direct access to the quantum component 
of the collective mode noise. 

In the weak tunneling regime, the switching of the de- 
tectors can be described in terms of the master equation. 
The cross-correlator, on the other hand, is found using 
the extended master equation formalism, which describes 
the evolution of the probability distribution in the mixed 
space: for the charge, transmitted through the QPCs, 
and the occupations of the two detectors' states. Cor- 
responding generalized transition matrix includes both 
the rates of the transitions of the detectors and the cu- 
mulant generators of the currents in the QPCs. The 
cross-correlator is then obtained from the largest eigen- 
value of the generalized transition matrix. We propose 
two methods of the evaluation of this eigenvalue, each 
method having its own advantages for different symme- 
tries of the transition matrix in the cases of the classical 
and quantum noise. 

To find the transition rates we calculate directly the 
time evolution of the detectors state to the fourth order 
of the perturbation theory. While, typically, this leads to 
well known divergences originating from the continuum 
of the environmental states, we avoid those divergences 
in a well controlled way by using the idea of the sepa- 
ration of time scales. We consider a time interval much 
shorter than a characteristic time of the evolution of the 
detectors, and much longer than the correlation time of 
the noise, that induces transitions. On this time interval, 
the solution of the master equation suggests a linear in 



time variation of the occupation probabilities, while the 
formal perturbation expansion of the transition proba- 
bilities, with respect to the tunneling Hamiltonian, gives 
rise to the quadratic in time terms. We observe that the 
quadratic in time terms are reducible, and thus cancel 
by taking the logarithm of the time dependent transition 
matrix. This leads to the final results for the elements of 
the transition matrix. 

By turn, perturbative calculation of the time evolu- 
tion of the transition probabilities requires averaging the 
products of two and four vertex operators over the fluc- 
tuations of the collective mode. In the limit of weak 
coupling of the detectors to the collective mode, the cu- 
mulant expansion has been used, and the cumulants are 
expressed through the components of the zero-frequency 
expansion of the current noise spectral functions. Such 
an expansion is justified in the limit of weak coupling, 
because the main contribution to time integrals for the 
elements of the transition matrix comes from long times. 
Note that the weak coupling resummed in this way, leads 
to the classical broadening of the quantum resonances at 
£j = in the transition rates of each detector separately. 
But due to the local nature of these resonances, they do 
not contribute to the cross-correlator. 

However, interaction of the detectors with the collec- 
tive mode leads to the effective non-local coupling be- 
tween the detectors. Already to the 2nd order in tun- 
neling, this effective coupling contributes to the cross- 
correlator a smooth background, which is induced by the 
classical noise, and thus does not represent interest. 

Finally, a more interesting effect appears in the 4th or- 
der of the perturbation theory with respect to tunneling 
Hamiltonian. Namely, at the degeneracy points ei = ±£2 
of the two detectors, a narrow non-local resonance arises. 
We show that when the effective coupling between two 
detectors is reduced, the main contribution to this non- 
local resonance arises from the quantum component of 
the noise source current. We estimate this effect, com- 
pare it to classical contribution to the cross-correlator 
and propose the range of parameters where the quantum 
contribution dominates. We show that the observation 
of quantum noise is feasible in the regime of weak cou- 
pling. Moreover, we argue that a sharp non-local reso- 
nance from the quantum noise may be well visible over 
the smooth background of the classical contribution, and 
has a characteristic energy dependence that can be used 
to distinguish it from the classical resonance. 
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